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1. Introduction 


Our purpose 

( 1 . 1 ) 


is to study the following system of partial differential equations: 

(i) ut + \o''^Uxx — fu-\-H{t,Ux,[mUx\) = 0<t<r, 0<x< 

(m) rrit — \cr‘^mxx — (G(t, Ux, [mux])m)^ = 0, 0<t<r, 0<x< 

{in) m{0,x) = mo{x),u{T,x) = ut{x), 0 < x < L 

(iv) u{t, 0) = m{t, 0) = 0, Ux{t, L) = 0, 0 <t <T 

(v) ^a^mx(t, L)-f G(t, Ux(t, L), lmux])'m(t, L) = 0, 0<t<T 


L 

L 


where T > 0 and L > 0 are given constants, mo and ut are known smooth functions, and H and G 
are defined below in Section o We mention for now that H and G depend on the variable mUx 
in a nonlocal way, in particular they are functions of 

L 

Ux{t, x)m{t, x)dx. 

System (jl.ip was introduced by Chan and Sircar in [9] to represent a mean held game in which 
producers compete to sell an exhaustible resource. Here we view the producers as a continuum 
of rational agents whose “density” is given by the function m{t, x) governed by a Fokker-Planck 
equation. Each of them must solve an optimal control problem corresponding to the Hamilton- 
Jacobi-Bellman equation (ll.l|) fii. Further details will be given below in Section fl.ll 


I 


Mean held games were introduced in mi EH] to describe differential games with large numbers of 
players represented by a continuum. Most recent results deal with models of the form 

Ut + \a‘^Uxx -ru + H{t, x,Ux) = V[m] 
mt - \(T‘^mxx - {G{t,x,Ux)m)^ = 0 

where V[m] is a monotone function. There have been a number of existence and uniqueness 
theorems proved when V[m] depends on m both locally [23l (H (H (H El EH ttS] and non-locally [6]. 


Both authors are grateful to be supported in this work by the National Science Foundation under NSF Grant 
DMS-1303775. 
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More recently some progress has been made toward similar results for cases where the Hamiltonian 
depends on m in a nonlinear way [laiiiKiT]. However, none of these results address systems where 
the coupling happens in the nonlocal part of the Fokker-Planck equation. 

Applications of mean field games to economics have attracted much recent interest; see mmm 
for surveys of the topic. The model we study here, which comes from [9], describes Bertrand or 
Cournot economic competition in the mean field limit (i.e. for a continuum of producers/consumers). 
It resembles the model proposed by Gueant, Lasry, and Lions to model oil production m, and it 
appears in a more complex and highly nonlinear form in m to describe the response of traditional 
oil producers to new technological developments such as renewable energy and “fracking”. In 
|I0j the authors point out, “There does not exist anything like general existence and uniqueness 
theorems for PDE systems of this kind.” This statement has inspired the present work, in which 
we prove existence and uniqueness for (HI]). 

We would like to mention in this context a recent result of Burger, et al. [3| which provides an 
existence and uniqueness theorem for a mean field games model of knowledge growth introduced by 
Lucas and Moll m- Their model also involves a coupled Boltzmann/Hamilton-Jacobi system of 
equations in which the coupling occurs through an integral over the space variable. This structure 
is natural for applications to economics, since aggregate quantities such as market price or total 
production are expressed mathematically as averages with respect to the density of agents. For 
this reason it is desirable to develop techniques to analyze PDE systems of this type. 

In this article we prove that under general conditions there exists a classical solution to dn]), which 
under a certain restriction is also unique. By “classical solution” we mean that the equations in 
dni) hold pointwise. We consider only the case where cj > 0 so that the equations are of parabolic 
type. Existence is obtained by applying the Leray-Schauder fixed point theorem; accordingly, the 
main effort of this paper is to provide a priori estimates of solutions. A new feature of our analysis 
is the estimation of the nonlocal term Ux{t,x)m{t,x). Although traditional methods provide 
estimates of u in L°°, it is not immediately clear how to obtain similar estimates for the gradient 
Ux■ In Section 12.31 we exploit the structure of (II.Hi by directly computing the time derivative of 
the nonlocal term, careful analysis of which allows us to derive higher order regularity. 

The remainder of this paper is organized as follows. In the rest of the introduction we give definitions 
of the functions H and G from CH), introduce some notation and give our main assumptions on 
the data. Section [2] is devoted to a priori estimates and constitutes the core of this paper. In 
Section [3] we prove the existence of solutions. Finally, in Section H we prove uniqueness under an 
additional hypothesis. 

1.1. Specification and explanation of the model. We summarize the interpretation of (|l.ip 
as follows. Let t be time and x be the producer’s capacity. We assume there is a large set of 
producers and represent it as a continuum. We say m{t, x) is the “density” of producers at time f, 
so that 

(1.2) r]{t) := f m{t,x)dx, 0 < r]{t) < 1 

Jo 

represents the total mass of producers remaining with positive stock. Note that r]{t) is a decreasing 
function in time. 

The first equation in (11.11) is the Hamilton-Jacobi-Bellman (HJB) equation for the maximization 
of profit. Each producer’s capacity is driven by a stochastic differential equation 

(1.3) dX{s) = -q{s,X{s))ds + alx{s)>odW{s), 

2 


where q is determined by the price p through a linear demand schedule 
(1.4) q = (p, p) = a{r]) — p + c{p)p, p > 0. 


In (|1.4I) p represents the market price, that is, the average price offered by all producers. This is 
given by 


(1.5) 




where p*{t,x) is the Nash equilibrium price. The coefficients 0 ( 77 ) and c(r/) = 1 — 0 ( 77 ) are defined 
by 


( 1 . 6 ) 


0 ( 77 ) := 


1 

1 + 677’ 


0(77) : = 


67 / 

1 + 677 


for a given fixed competition parameter 6 > 0. The case 6 = 0 corresponds to monopoly, while 
perfect competition is given by 6 = + 00 . Thus each producer competes with all the others by 
responding to the market price. 


We define the value function 

(1.7) u{t,x) :=supe|^ e~'^^^~^"^p{s)q{s)ds + ut{X{T)) \ X{t) = x 

where q{s) is given in terms of p{s) by (II.4p . The optimization problem (II.7p has the corresponding 
Hamilton-Jacobi-Bellman equation 

(1.8) Ut + Ircr'^Uxx -ru + max [{a{p{t)) -p + c{p{t))p(t)) (p - u^)] = 0. 

i P 

The optimal satisfies the first order condition 

(1.9) p*{t, x) = ]^ {a{p{t)) + c{p{t))p{t) + Ux{t, x)) , 
and we take q*{t,x) to be the corresponding demand 

( 1 . 10 ) q*{t, x) = ]^ (aivit)) + c{p{t))p{t) - u^it, x)). 

This leads to Equation (jl.ip fil where 

(1.11) H{t, Ux, [mUx]) := q*{t, xf = ^ ia{p{t)) + c{p{t))p{t) - Uxf 

On the other hand, the density of producers is driven by the Fokker-Planck equation (jl.ip (ii). where 

( 1 . 12 ) G{t, Ux, [mUx]) := q*{t, x) = ^ {aivit)) + c{p{t))p{t) - Ux) 

The coupling takes place through the average price function, which, thanks to (II.5p and (jl.Op . is 
given by 

(1.13) p{t) = ^ x)m{t, x)dx^ 

We have taken Dirichlet boundary conditions at x = 0 as in [9]. On the other hand, rather 
than taking L = +00 and working on an unbounded domain, we have taken Neumann boundary 
conditions at x = L, which represents a diffusion which is reflected at this boundary point. We can 
think of L as an upper limit on the capacity of any given producer. 
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1 . 2 . Notation and assumptions. Throughout this article we define Qt ■= (0,7") x ( 0 ,L) to be 
the domain, St ■= ([0,7"] x {0,L}) U ({T} x [0,L]) to be the parabolic boundary, and at times 
Tj’ := ([0,T] X {0}) U ({T} x [0, Lj) to be the parabolic half-boundary. For any domain X in M or 
we dehne LP{X), p £ [1, -|-oo] to be the Lebesgue space of p-integrable functions on X; C^{X) 
to be the space of all continuous functions on X; C'“(X), 0 < a < 1 to be the space of all Holder 
continuous functions with exponent a on X; and C^~^°‘{X) to be the set of all functions whose 
n derivatives are all in (^"(X). For a subset X C Qt we also define C^’^(X) to be the set of 
all functions on X which are locally continuously differentiable in t and twice locally continuously 
differentiable in x. By (^“/^’“(X) we denote the set of all functions which are locally Holder 
continuous in time with exponent a/2 and in space with exponent a. 

We will denote by C a generic constant, which depends only on the data (namely ut, m-o, L,T,a,r 
and e). Its precise value may change from line to line. 

Throughout we take the following assumptions on the data: 

(1) ut{x) and mo{x) are functions in C'^''''>'([0, Lj) for some 7 > 0. 

(2) ut and mo satisfy compatible boundary conditions: ut(0) = u'j,{L) = 0 and mo(0) = 
mo(T) = mo(L) = 0. 

(3) mo > 0 and moix)dx = 1, i.e. mo is a probability density. 

(4) Ut >0 and u'rp > 0, i.e. ut is non-negative and non-decreasing. 

Remark 1.1. Of all the assumptions, the stipulation that ut be non-negative and non-decreasing 
seems the least essential; it is not necessary for most estimates. However, it appears to he needed 
to prove the a priori bounds of Section \2.3[ 


2. A PRIORI ESTIMATES 

The goal of this section is to estimate various norms of solutions to using constants depending 
only on the data. In Section 12.11 we prove some standard results, including the usual “energy” 
type estimate on the quantity Jq /q Ux''^ dxdt. Then in Section \2^ we prove a priori bounds on 
the solution to the Hamilton-Jacobi equation using classical techniques for parabolic equations. 
Section o is our most original contribution; there we show that the term jJ" Uxit, x)m{t, x)dx is a 
priori bounded uniformly in t. Finally, in Section [2.41 we use the previous estimates to prove higher 
regularity. 


2.1. Basic a priori estimates. 

Proposition 2.1 (Main a priori estimates). Suppose {u, m) is a pair of smooth functions satisfying 
(|1.1[) . Then 

(2.1) u{t,x) > 0,m{t,x) > 0, ||m(t)||ii(o,L) < ||mo||i,i(o,L) Vt G [0,r],Vx G [0,L]. 

Moreover, for some C > 0 depending on the data, we have 

I / mu^dxdt < C, 

0 Jo 



m\G{t,Ux, [mux])\‘^dxdt 



mH{t,Ux, [mux])dxdt < C. 


( 2 . 2 ) 

and 

(2.3) 
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Proof. We start by proving (12.11) . Let 4> be any function in C'^(M) with = 0. Multiply (ll.lM iil 
by and integrate by parts to get 


(2.4) f 4>{m){t,x)dx = --a'^ f f cj)"{m)m‘^dxds + ( f 4>"{m)mxGm-dxds. 

Jo ^ Jo Jo Jo Jo 

Take 4>{s) = where s_ := (|s| — s)/2 and let 5 —)• 0 to deduce 

nL rt pL pt pL 

(2.5) / m_(t,x)^(ix = —/ / \{m-)x\^dxds+ 2 / / {m-)xGm-dxds 

Jo Jo Jo Jo Jo 


< - 


[ [ |(m_)^pdxds +[ [ {m-fdxds, 

2 Jo Jo Jo Jo 


where we note that (mo)_ = 0. By Gronwall’s inequality we obtain m-{t,x) = 0, which proves 
positivity. On the other hand, if we take 4>{s) = and let 5 —)• 0, then we deduce ||?Ti(t)||Li(o,i,) < 

ll^-ollLqo.L)- 

Also, since we have 

( 2 . 6 ) 


a 


-ut- —Uxx + ru > 0, 

we can deduce using similar arguments that u > min^, tt(T, x). Hence, in particular, we have 
u > 0 by the assumption ut > 0- Thus we have proved (12.ip . 

Next we prove (j2.2p . from which follows ()2.3jl . Multiply (jl.ip (il by m and (jl.ip iii) by u and 
integrate by parts to get 


pLi pLi pi pL/ 

(2.7) / UT{x)m{T,x) dx — / u{0,x)mQ{x) dx = r / / um dxdt 

Jo Jo Jo Jo 

rT pL pT rL 

— / mH{t,Ux,[mUx]) dxdt — / / 'mUxG(t,Ux,[mUx]) dxdt. 

Jo Jo Jo Jo 


Since u > 0 and m > 0, we get 

( 2 . 8 ) 


pL/ pT pL 

/ UT{x)m{T,x) dx + / / m[uxG{t,Ux,[mUx]) + H{t,Ux,[mUx])] dxdt > 0, 

Jo Jo Jo 

and then using the fact that ||m(t)||^i < 1 we can rewrite this as 

pT pL pT pL pT 

(2.9) / / mu^ dxdt < / / m(a + cp)^ dxdt + 4||ttT||oo = / ri{t){a + cp)‘^dt + 4:\\ut\\ 

Jo Jo Jo Jo Jo 

since a + cp does not depend on x. 

To analyze the right-hand side, we observe that 

2 


a{r]{t)) + c{r]{t))p{t) = 
By Cauchy-Schwartz we see that 




2 eri{t) 2 + er/(t) Jo 
e?7^/^(t) 


Ux{t, x)m(t, x)dx. 


\a{r]{t)) + c{r]{t))p{t)\ < 1 + ul{t,x)m{t,x)dx'^ 
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1/2 


( 2 . 10 ) 







which implies, using the fact that 0 < ? 7 (t) < 1, 


r]{t){a{r]{t)) + c{r]{t))p{t))‘^ < (1 + + (1 + ul(t, x)m{t, x)dx 


< 1 + 1/5 + 


(l + 5)e 


2 + e 

for an arbitrary 5 > 0. By choosing 5 = 1/e, then (j2.9|) becomes 

r-T rL 


/ u^{t,x)m{t,x)da 

Jo 


f f dxdt < (2 + e)(l + e)T + 4(2 + e)||riT| 
Jo Jo 


which yields (j2.2p . As for (j2.3p . we combine fj2.2p with (|2.lOp and the definition of G and H. EH 


We may now deduce certain a priori bounds on the Fokker-Planck equation, which will be useful 
later on. 

Lemma 2.2 (Regularity of m). Suppose {u,m) is a pair of smooth functions satisfying (|l.lh . 
Then there exists a constant C > 0 depending on the data such that 


rT rL 


( 2 . 11 ) 


'0 


/ 


mz 


m + 1 


-dxdt < C. 


Proof. Multiply (ll.lU fiil by ln(m + 1) and integrate by parts to get 

rL ^rL 

2 Jq 1 + m 


(2.12) ^ j (j)im){t) dx 


Gmrrix . 

^ -dx - / - -dx 


/o 1 + m 


m't 1 f G^m^ 

A Jq 1 + m Jq 1 + m 


L 


< - 


where (f{m) = (1 + m) ln(l + m) — m. By Equation 12.31 in Proposition 12.11 and using the fact that 
mo is bounded and > 0, we get 


(2.13) 


,2 rT rL z rr i rl rl^ 

/ / - —dx < / 0(mo) dx -\—^ / G^m dxdt < G. 

/o /o 1 + ^ Jo Jq Jq 


2 rL 

A^dx < 

+ m 


T rL 


□ 


2.2. A priori bounds for the Hamilton-Jacobi equation. Let f{t) := a{p{t)) + c{rj{t))p{t). 
Then (ll.lD lil reads as 

(2.14) - rtt + ^(/(t) - = 0, 

from which we can estimate 

(2.15) -ut-+ ^/(t)^ + iu^. 

From Proposition 12.II we know that / € L^(0, T) with an a priori bound on its norm. Using classical 
arguments, this is enough to infer an estimate on u as well as an estimate on Ux, as the 
following proposition makes clear. 
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Lemma 2.3. Suppose u is a smooth function on [0, T] x [0, L] satisfying 

(2.16) -ut- kuxx < g{t) + jul, u{T, x) = UT 

where j,k are positive constants, g{t) > Q is an integrable function on [0,T], and ut = ut{x) is 
a smooth function on [0,L]. Assume that u is bounded below, that u{t,0) = 0, and Ux{t,L) = 0. 
Then there exists a constant C = C{j,k,UT, ||ff||Li(o,T)) such that 

(2.17) Ikiloo + [ [ Ux dxdt < C. 

Jo Jo 


Proof. Let w{t, x) = exp |a (^{t, x) + Jq g{s)d^ | — 1 for A = j'/fe. Then 

(2.18) — wt — kwxx < (j — kX)Xe^'^Ux = 0. 

In particular w satisfies the maximum principle, i.e. 

max w = maxtc 
[0,r]x[0,L] Ft 

where = ([0, T] x {0})U ({T} x [0, L]). To see this, it suffices to take // = maxr^, w, then multiply 
(j2.18p by {w — ia)+ and integrate by parts (note that Wx{t,L) = 0) to get 

rL pT rL rL 

/ {w{t,x) — ii)\dx = —k / / {{w — p)j,.)‘l,dxdt-\- / {w{T,x) — p)\dx <t). 

Jo Jo Jo Jo 

It follows that w < fi everywhere, from which we deduce that < y ^In the other hand, by 

the definition of w we can directly compute 

rT 


/i = max exp 

x^[0,L] 


X ^t{x) + J g{t)dt 


-L 


which is a constant depending only on Hurlloo and 11 ( 7111,1 (o,r)- 
Using the same computation with w instead of (rc — //)+ we get 


ff 


wldxdt < 


irj, ^2 , ^ 

w[l , x),dx < ——. 
1 k 


Since Wx{t,x) = Aexp |a (^u{t,x) + /q 5 '(s)(is^ | Ux{t,x) and u is bounded below, we deduce 


ff 


Ux dxdt < 

“ X^k 


□ 


2.3. Analysis of nonlocal term. In order to obtain higher regularity on u, we need to analyze 
the nonlocal coupling term Ux{t,x)m{t,x)dx. In particular, we will show it is bounded. 


In the case when a = 0, we have a fortuitous identity which follows from integration by parts. 
Differentiate (jl.lip to get 

Uxt -rux-^{a + cp- Ux)uxx = 0 , 

noting that a = 0. Then multiply by m and ()1.12l) to get 


jL 

dt 


f 


Ux{t, x)m{t, x)dx = r 
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Ux{t, x)m{t, x)dx, 










which means 



Thus, as long as ut is smooth, we know that the term Ux{t,x)m{t,x)dx is bounded uniformly 
in t. This in turn implies that p{t) is bounded, which allows us to analyze the regularity of u by 
classical methods for parabolic equations. 


Unfortunately, when u > 0 this formal calculation fails; we get instead 



0 


with no estimates on the boundary terms. On the other hand, thanks to the following lemma, we 


note that each of the terms Ux{t, 0)mx{t, 0) and Uxx{t, L)Tn{t, L) has a definite sign. This will allow 
us to prove that each of these terms is integrable in time with an a priori bound on its L^(0,T) 
norm. 

Lemma 2.4. Suppose (u, m) is a pair of smooth functions satisfying dni). Then Ux{t,x) > 0, 
"rx(LO) > 0, and Uxx{t,L) < 0 for all {t,x) G [0,T] x [0, L]. 

Proof. Notice that, by Proposition 12.11 u and m are both non-negative, hence their minimum is 
attained at u{t,0) = 0 and m(t,0) = 0, respectively. It follows that Uxit,0) and mxit,0) are 
non-negative for all t £ [0, T]. 

Differentiate ow to get 

(T^ 1 _ 

(2.20) Uxt H ~Uxxx '^Ux 2 ^xjUxx — 0- 

Note that Ux{t, 0) and 0 = Ux{t, L) are both non-negative. As u'rp{x) is also non-negative, it follows 
that Ux > 0 everywhere by a classical maximum principle argument. Thus Ux{t,L) = 0 is a 
minimum of Ux, so Uxx{t^L) <0. □ 

Returning to (j2.19p . we note that the terms on the right-hand side have definite but opposite signs. 
To show that each of them is integrable in time, we localize away from each boundary point. Once 
this is accomplished, we can see that the term Ux{t,x)m{t,x)dx is bounded. We prove this in 
the following: 

Proposition 2.5. Suppose {u,m) is a smooth solution of O)- 
(i) For any 5 > 0, there exists a constant Cs > 0 such that 



■L-5 


( 2 . 21 ) 



(a) For any 5 > 0, there exists a constant Cs > 0 such that 



( 2 . 22 ) 


(Hi) By (i) and (ii), there exists a constant C depending only on the data such that 



(2.23) 
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Proof. Take a smooth, non-negative function ^ = ^(x) on [0, L] to be further specihed later. Mul¬ 
tiply (12.201) by and integrate by parts using (ll.llh iil to get 


(2.24) e'’*^ J u^{t)m{t)C, dx^ = -^CWu^rxit, L)m(t, L) - YC(0)ux(t,0)m:r(t,0) 


CxUxmx dx — 


^2 rL 1 1 

/ CxxUxm dx --{a + cp) / (xUxm dx + - / Cxulm dx. 

^ Jo ^ Jo ^ Jo 


Let us estimate the time integral of each of the terms in the second line of (j2.24p . First we have 

' y J ^ ^ ' y i ^ ^ ^ 2 •v 

(2.25) i rf Cx^x'^x dx \ dt ^ <7 11 Coil I oo 11 { ul{m + 1) + I dxdt < CIlCxIlcx)- 
using Proposition 12.II and Lemma 12.31 Likewise, 


(2.26) 

(2.27) [ 

Jo 

and finally 

(2.28) 


,2 rL 


- [ 
2 Jo 


CxxUxITl dx dx 


pT rL 


{a + cp) / CxUxTLl dx 

Jo 


dt< IICxxllcxD / / {(^x T l)m} dxdt < CIlCxxllc 

^ Jo Jo 

hi < IlCxIloo j ^{a + cp)^ + j dxj df < CIICx 


/ / CxU^m dx 

Jo Jo 

To summarize, we may write (12.241) as 
(2.29) 


dt < IlCsIloo / / ulm dxdt < C\\Cx 

Jo Jo 


^rtfj_ ( e J Ux{t)m{t)C dx^ = -^C,{L)uxx{ti L)m{t, L) - ^C,{t))ux{t, t))mx{t, 0) -b If (t), 


dt 


where 


r\ic{m<c{c) 

Jo 


such that C'(C) depends only on ||Cx||oo) ||Ca:a:||cx)) and previous estimates. 

Now let us prove (i). We specify that C,{x) = 0 for L — 5/2 < x < L, ({x) = 1 for 0 < x < L — 5, 
and 0 < C < 1- Then we can assume C{C) < Cs, where Cs is some constant proportional to 
for d > 0 small. Integrating (j2.29p over [0, T] we get 


a 


2 rT 


(2.30) — e ^^Uxit,0)mxit,0)dt 

2 Jo 


rL rL ri 

= / Ux{0,x)mo{x)C dx — e~^'^ / Uj^{x)m{T,x)C dx + / e~^^I(^{t)dt. 

Jo Jo Jo 

Now on the one hand, using the fact that Ux > 0 and that u is bounded lLemma l2.3p . we have 
(2.31) 

rL rL rL 

/ |ttx(0,x)mo(x)C| dx = / Ux{0,x)mo{x)C dx < \\mo\\oo na;(0,x)dx = ||mo||ooii(0,L) < C. 
do Jo Jo 

On the other hand, since fj' m(t, x)dx < 1 for all t, we have 

/ |u^(x)m(r,x)C| dx < llu^lloo- 
do 


(2.32) 



















Using the fact that Ux >0 and mx{t,0) > 0 from Lemma 12.41 we deduce that 

l-T pT 

(2.33) / \ux{t,0)mx{t,0)\dt < / e^^Ux{t,0)mx{t,0)dt < Cs- 

Jo Jo 

Finally, integrating (|2.29l) this time over [t,T] we get 

(2.34) / Ux{t,x)m{t,x)C dx 

Jo 


2 rT 


pi pL, pi 

/ e~^^Ux{t,0)mx{t,0)dt + e~^'^ / u'rp{x)m{T,x)( dx — / e~^^Ii^{t)dt 

Jt Jo Jt 

from which we obtain (j2.21jl . 

In a similar way we can prove (ii). This time we specify that C(x) = 0 for 0 < x < 5/2, Ci^) = 1 
for 5 < X < L, and 0 < C < 1- Again we can assume C{(/) < Cs, where Cs is some constant 
proportional to 1/5^ for <5 > 0 small. Integrating (I2.29P over [0, T] we get 


(2.35) - 


cr 


2 rT 


e ^^Uxx{t,L)m{t,L)dt 


/ L pL pT 

Ux{0,x)mo{x)C dx + j u'rp{x)m{T,x)Q dx — j e~^^I(^{t)dt. 


JO JO 

Now since Uxxit,L) < 0 from Lemma 12.41 and m > 0, we can deduce 

(2.36) [ \uxxit,L)m{t,L)\dt <Cs 

Jo 

using the same estimates as in the proof of part (i). Now integrating (12.291) over [t,T] we get 

(2.37) f Ux{t,x)m{t,x)( dx 

Jo 


a 

'2 


2 rT rL rT 

e~'^^Uxxit, L)m{t, L)dt + / u'rp{x)m{T, x)</dx — / e~^^lQ{t)dt 

Jo Jt 


from which we infer (|2.22p . 

Finally, p2.23p follows from p2.2ip and (12.221) by fixing any 5 < L/2. □ 

Corollary 2.6. Suppose {u, m) is a smooth solution of ([HI). Then for some constant C depending 
only on the data, 

(2.38) Hdit)) + c{vit))Pit)\ <C yt€ [0,T]. 

Proof. By the definition of a, c, and p, it is enough to have \ Ux{t,x)m{t,x)dx\ < C for all 
t G [0,r], which we get from Proposition 12.51 □ 

2.4. Full regularity of u. Let us return to ()2.14p : 

ut + '^Uxx -ru + ^ifit) - Uxf = 0, 
where we recall f{t) := a{r]{t)) + c{rj{t))p{t) . We can now write 


(2.39) 


a 


1 


-Ut - —Uxx + ru<Ci + 
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where Ci is the constant coming from Corollary 12.61 It is now possible to obtain global estimates 
on \ux\, which we will then be able to “bootstrap” to gain higher regularity on u. 

Lemma 2.7. For {u,m) a smooth solution of (ini), we have 

(2.40) \ux{t,x)\ < C V(t,x) G [0,T] X [0,L]. 

for some constant C depending only on the data. 


Proof By Lemma 12.41 we already know Ux > 0. It suffices to prove that Ux < C. We claim 

(2.41) Ux{t,0)<C. 

To see this, set v = e^/^^ — 1 and use (j2.39l) to get 

(2.42) -vt- ^Vxx < =: A- 

2 

Note that v{t, 0) = Vx{t, L) = 0. Set v = v + Me~^ where M is large enough that 

Vx{T,x) = \e^-r^^^/^\'T{x) < Me-^ 

(7^ 

and 

2 

Cl < 

Then, on the one hand, we have —vt — ^Vxx < 0, and since Vx{t, L) = —Me~^ < 0 we get as before 
the maximum principle 

max V < maxu. 

[0,T]x[0,L] Ft 

On the other hand, we have Vx{T, x) < 0 and so v{T, x) < v{T, 0) = M for all x G [0, L], It follows 
that v{t,0) = M is the global maximum of v, hence Vx{t,0) < 0 at each t G [0,T]. Recalling the 
dehnition of v we get 

-^e^/'^^Uxit,0) < M, 

and since u > 0 we get Ux{t, 0) < Mcr^, which is the desired estimate. 


Taking into account the assumption that ut is smooth, we have thus shown 
(2.43) max|u 3 ;| < C, Ty ;= ([0,T] x {0,L}) U ({T} x [0,L]) 


where C depends only on the data. 

Now differentiate (j2.14p to get 

(2.44) Uxt "^Uxxx T FUa; ~{f Ux')Uxx — 0- 
Then w{t,x) = Ux{t,x)e~^^ satisfies 

(2.45) - wt - + \if ~'^x)wx = 0- 

By classical arguments, w satisfies the maximum principle, i.e. 

\w{t,x)\ <max|r(;| <max|ua;|, 
r 7^ r 7^ 


from which it follows that 


Ux{t,x)\ < e”'^max|ua;| < C. 
Ft 


□ 
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Corollary 12.71 permits us to obtain higher order estimates for u. 

Proposition 2.8. There exists a constant C depending only on the data such that if {u, m) is a 
smooth solution of (ini), then for some a > 0 

(2.46) \\u\\(jl + a/2,2+OL((^^ + \\'m\\(jl + a/2,2 + a (^) ^ C 

where is the parabolic Holder space defined on Qt = [0,T] x [0, L]. 

Proof. Observe that u and Ux each satisfy a linear parabolic equation with coefficients which are 
bounded by constants depending on the data. By [201 Theorem IV.9.1], we have that u is bounded 
in LP{0, T; W'^’P{0, L)) CiW^’P(Qt) for any p > 1 and thus u, Ux and Uxx are all bounded in a Holder 
space C^{Qt) for some (3 > 0. Then by |20l Theorem III.7.1] m has an a priori bound in L°°{Qt). 
Further, we observe that (jl.lh fiil can be written 

1 _ 1 

(2.47) mt - —mxx - i^[.a + cp - UxjrUx + -Uxxm = 0, 

which also has coefficients bounded by the data. Using the same technique as in Lemma 12.71 we 
obtain an a priori estimate on mx{0,x). Then (|2.19l) can be used to directly estimate the term 
Jq Ux{t,x)m{t,x)dx in C'^([0,T]). We can now see that (ll.llh il.fiil are both parabolic equations 
with coefficients estimated in Holder spaces by constants depending only on the data. Applying 
|20l Theorems IV.5.2 and IV.5.3] now gives the conclusion. □ 


3. Existence 


We now prove the main result of this paper. 

Theorem 3.1. There exists a classical solution of (11.11) . 

Proof. We use the Leray-Schauder hxed point theorem. Consider the map {u, m) e->■ (u, /) = 
T{u,m-,T) given by solving the following parametrized set of PDE systems: 



o<t<r, o<x<L 
0 < t < T,0 < X < L 
0 < X < L 
0<t<T 
0<t<T 


(3.1) {in) f{0,x) = Tmo{x),v{T,x)=TUT{x), 


{iv) v{t, 0) = f{t, 0) = 0, Vx{t, L) = 0, 

{v) \a‘^fx{t, L) + rG(t, Vx{t, L), [mvx])f{t, L) = 0, 


for r G [0,1]. Let X be the space of all {u,m) such that u and Ux are both Holder continuous. 


say in C^{Qt), and m is in IU^’°°(0, T; L^(0, L)). Then by inspecting the dehnitions of G (|1.12j) 
and H (jl.llll we find that H{t,Ux, [mUx]) and G{t,Ux, [mUx]) are both Holder continuous as well. 
By m Theorem IV.5.2] there is a solution v of (|3.ip (i) satisfying (hi) and (iv) such that v G 
(ji+a/ 2 , 2 +a^Q^'^ fgj. gome o > 0 and 

(3-2) llull(^i+c,/2,2+c,(^) < G\\{u, m)\\x = G{\\u\\^a/2,c(Q:^) + llwa;llc«/2.c(^) + ll"i|llTl.°°(0,r;Ll(0,L))) 
for some generic constant G. Next, we write (I3.1|) ('iil as 



and we note that the coefficients are Holder continuous. Furthermore, because m G IU^’°°(0, T ; L^(0, L)) 
we can see that 
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is independent of x and has a bounded time derivative. Then we can apply [201 Theorem IV.5.3] 
to get a solution / e satisfying ()3.1l) fiiil and (v) such that 

(3-3) ||/||cl+a/2,2+c«(^) < C'(||t’||cl+c./2,2+a(^) + 11 11 (0,T;Ll (0-.^-)) - ^ I11 ^' 

It follows that T : V X [0,1] —)■ X is a well-defined and compact mapping, since x 

(ji+a/ 2 , 2 +a jg compact in X by the Arzela-Ascoli Theorem. 


Now for r = 0 we have T(u,m-,0) = 0 for all {u,m) by standard theory for the linear heat 
equation. On the other hand, if {u,m) G X, r £ [0,1] is such that {u,m) = T{u,m-,T) then 
{u,m) G X is a solution of (11.11) with mo,UT, G and H replaced 

by TmQ,TUT-, tG and tH, respectively. Then the a priori estimates of Section [2] carry through 
uniformly in r G [0,1], and so by Proposition 12.81 we obtain a constant Cq depending only on the 
data such that 


\\{u,m)\\x < Gq. 

By the Leray-Schauder fixed point theorem (see e.g. m Theorem 10.6]), there exists {u,m) such 
that {u,m) = T{{u,m)‘, 1), which is a classical solution to (II.ip . □ 


4. Uniqueness 

The structure of (jl.ip makes uniqueness a nontrivial issue. Unlike traditional mean field games in 
which uniqueness is verified by a straightforward use of the “energy identity” (thanks to the fact 
that the coupling is monotone EH), System (11.11) does not allow such an argument. Our uniqueness 
result will rely heavily on the fact that solutions are smooth with a priori bounds, and it will hold 
only when e is small. 

We now proceed to state and prove the main uniqueness result. 

Theorem 4.1. There exists eo > 0 sufficiently small such that for any e < cq , (HI) has at most 
one classical solution. 


Proof. Let and {u 2 ,m 2 ) be two solutions. Define, for i = 1,2, 

1 1 

Hi = -(a(? 7 j(t)) + c{r]i{t))pi{t) - dxUi{t,x))'^, Gi = -{a{r]i{t)) + c{rii{t))pi{t) - dxUi{t,x)) 

where pi{f) and pi{t) are defined according to the definitions in (II.2p and (jl.l3p . with u and m 
replaced by Ui and m,. Then, in particular, u = ui — U 2 satisfies 

2 

(4.1) Ut + ^Uxx - ru + Hi - H 2 = 0, u{T, •) = 0 

while m = mi — m 2 satisfies 


(4.2) 


cr 


mt - —mxx - {Gimi - G2m2)x = 0, m(0, •) = 0. 


Let us introduce some notation. Observe that 

2 


aiViit)) + c(r?i)(t))pi(t) = 


+ 


2 + epi{t) 2 + €pi{t)jQ 


/ dxUi{t,x)mi{t,x)dx. 

Jo 


With this in mind, we define 


2 e f 

= TT, -TTT, Bi{t) = — - — dxUi{t,x)mi{t,x)dx. 

2 + epi[t) Jo 


2 + epi{t)' 


Notice that 2Gi = Ai + Hi — dxUi and Hi = Gf. 
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Now multiply (I4.ip by mi — m 2 and (|4.2p by ui — U 2 , integrate by parts and add to get the typical 
energy identity for mean field games: 

(4.3) f f e~^*[{Hi - H2){mi - m 2 ) + (Gimi - G2m2){dxUi - dxU2)]dxdt = 0, 

Jo Jo 

which can be rearranged to get 

rT rL 


(4.4) f f ^ '^*mi[Hi — H 2 + {dxUi — dxU2)Gi]dxdt 
Jo Jo 

+ [ [ e~^*m2[H2 - Hi + {dxU2 - dxUi)G2]dxdt = 0 . 

Jo Jo 


This, in turn, can be rearranged to give 

rT 


(4.5) -f e-^^{Ai{t)^-A2{t)^){vi{t)-m{t))dt 
Jo 

+ / e~'"^{mi + m2){dxUi — dxU2)‘^dxdt 

Jo Jo 


= h+l2 


where 


(4.6) h := [ e-"*(r/i(t) - m{t)){Bi{t)^ - B2{t)^ + 2Ai{t)Bi{t) - 2A2{t)B2{t)}dt 

Jo 

and 

(4.7) 


/■T rL 

h :=2 

'0 Jo 

By a simple computation we have 

l-T 

(4.8) 


[ [ e '^^{m 2 dxUi - midxU 2 ){Ai{t) + Bi{t) - A 2 {t) — B 2 {t))dxdt. 

Jo Jo 


8e 


^ e (^i(t)^ - A 2 {t)'^) imit) - m{t))dt > ^"^^3 / e ''\T]i{t) - Tj 2 {t))'^dt 


so we can write 

rT 


(4.9) 8e [ e - r] 2 {t))‘^dt + [ [ e "'^{mi + m 2 ){dxUi - dxU 2 )‘^dxdt 

Jo Jo Jo 

Our main task is to estimate Ii and I 2 ■ For this we will use the a priori estimates from Section [2l 
which say in particular that \dxUi\ < Cq for some Gq depending only on the data and on cq, where 
e < cq. Since we are going to make e small, we need to keep in mind that the constant Gq does not 
change for decreasing values of e. 

First, we estimate \Bi{t) — B 2 {t)\. We have 

(4.10) Bi{t)-B2{t) = 

■ / dxUi{t,x)mi{t,x)dx - ^ / dxU 2 {t,x)m 2 {t,x)dx 

Jo 2 + 6772 (t) Jo 


< Jl + 12- 


2 + 67?! (t) 


+ 


2 + 6772 (t) 

L 

dxUi (t, x) 77 ii (t, x)dx 


2 + e77i(t) 2 +er] 2 {t) J Jii 

fL 

dxiui{t, x) — U 2 {t, x)){mi{t, x) + m 2 {t, x))dx 


2(2 + €772 (t)) Jo 


+ 


2(2 + er]2{t)) Jo 


/ {dxUi{t,x) + dxU 2 {t, x)){mi{t, x) - m 2 {t,x))dx. 
Jo 
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Observe that 


2 + er]i{t) 2 + er] 2 {t) 


< 


e\r]i{t) - mit)\ 


j \mi{t,x) — m 2 {t,x)\dx < x) — m 2 (t,, and 

J (dxUi + dxU2){mi - m2)dx, < \f2 {d^ui + dxU2)‘^imi - m2)da 


1/2 


by the Cauchy-Schwartz inequality. So using the uniform pointwise bounds on Ux{t,x), we obtain 

1/2 

{dxUi - dxU 2 )‘^{mi + m 2 ) dx ] 

/o 


(4.11) \B^{t)-B2{t)\< 


2y/2 


{dxUi - dxU 2 )‘^{mi + m 2 ) dx 

+ - ^ 2 (t)| + \mi - 


1/2 


On the other hand, we have \Bi{t)\ < Coe/2 for z = 1,2. We deduce from (14.lip that 

1/2 


(4.12) \Bi{t)'^ - B 2 {t)‘^\ <-^^Co {dxUi-dxU 2 )‘^{mi+m 2 ) dx 

+ ^Colmit) - ^ 2(01 + \mi - m 2 |^dx^ 


1/2 


\Ai{t) - A2{t)\ < -\rj2it) - Vi{t)\- 


We also have 

(4.13) 

as well as |2lj(t)| < 1 for i = 1, 2. Combine this with (|4.1ip to get 

(4.14) \Ai{t)Bi{t) - A 2 {t)B 2 {t)\ <{dxUi - dxU 2 )^ {mi + m 2 ) dx 


1/2 


+ j(l + Co)|?7i(t) - r]2{t)\ + |"ii - m 2 \‘^dx 

From ()4.12p and (j4.14p it follows that 


1/2 


(4.15) Ii< 


e^Co + e r 


2y/2 




1/2 


e iviit) - mit)) / {dxUi - dxU2) (mi + m 2 ) dx dt 




+ 


e^Co + €^(1 + Co) 


c^Cq + eCo 
4 

rT rL 


f e ^\r]i{t) - r] 2 it))‘^dt+ 

Jo 

L^/^y e"^‘(7/i(t) - r/2(t)) ^y |mi-m2|^iix^ 


1/2 


dt 


pi pL/ pT pL 

(e) / / e“'’*(mi + m2)(i9a;tti — cl,rn 2 )^(ix(it + P 2 (e) / / e~''^{m 2 — mi)"^dxdt 

Jo Jo Jo Jo 


<Pi 
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where Pi(e),P 2 (e) —)• 0 as e —>■ 0. As for / 2 , setting D{t) = Ai{t) + Bi{t) — A2{t) — B2{t), we write 

(4.16) I 2 = [ [ D{t){m2 — mi){dxUi + dxU2)dxdt 

Jo Jo 

+ / / e~^^D{t){mi + m2){dxUi — dxU2)dxdt 

Jo Jo 

rT / rL \ 1/2 

< 2 CoL^/^ e-"'|Zl(t)|n (m2 — mi)^dx\ dt 

{nil + m 2 ){dxUi - dxU 2 )‘^'^ dxdt. 


By (|4.1ip and ()4.13p we have 

fT rL 


rl rL/ rT rL 

(4.17) I 2 < Poie) / / e“^‘(m2 - mi)^dx(it + P4(e) / / e“'’*(mi + m2)(i9a:iii - f/xii2) 

./o ./o io Jo 


where Psie), P^ie) —)• 0 as e ^ 0. By (I4.15P and (I4.17p . Equation (14.91) becomes 

(4.18) f f e~^^{mi+m2){dxUi — dxU2)‘^dxdt<-^‘^^^'^~^^^^^^ 

Jo Jo 1 


-Pl(6)-P4(e) Jo 


[ [ e ^*(m2-mi)‘ 

Jo Jo 


‘dxdt 


dxdt. 


Now we consider the Fokker-Planck equation (14.2p . Recall that m = mi — m 2 . Multiply by m and 
integrate by parts to get 

(4.19) - J w?{t,x)dx = —— [ [ mx dxdt — f f mx{Gimi — G 2 ni 2 )dxdt 


to JO 
v2 ri rL 


0 JO 
t pL 


< - 


— [ [ mx dxdt H—^ / [ {Gimi — G 2 m 2 )^dxdt 

4 Jo Jo Jo Jo 


t rL 


t rL 




|Gi — G 2 fm^ dxdt 


Recall that Gj = 4(Aj + Pj — dxUi). Then, using (j4.1ip . (j4.13p . and the fact that Gi and m 2 are 
bounded by some constant depending on the data, we obtain 


rL/ rt rL rt rL/ 

(4.20) / {mi{t,x)—m 2 {t,x))‘^dx<G / / {mi—m 2 )^dxdt+P^{e) / / {dxUi—dxU 2 )^m 2 dxdt 

Jo Jo Jo Jo Jo 


t pL 

0 Jo 


where P 5 (e) —)> 0 as e —)• 0. By Gronwall’s Lemma, we have 

rL rT rL 

(4.21) sup / {mi{t,x) — m 2 {t,x))‘^dx < P^{e) / / {dxUi — dxU 2 )‘^m 2 dxdt, 

te[o,r] Jo Jo Jo 

and then by appealing to (I4.18jl we have 

e-n(e)(P40 + P3(0) 

Jo Jo 


T 1 

(4.22) sup / {mi{t,x) — m 2 {t,x))'^dx < -- 

iG[0,T] Jo 1 


-Pi(e)-P4(e) Jo 


dxdt 


^ rj^^n(£M(£)±L(£)) sup 


1 - Pi(e) - P4(e) t€[o,T]. 


Fix e small enough; then (|4.22l) implies 


(4.23) 


sup / {mi{t,x) — m 2 {t,x))'^dx = 0, 
t£[0,T] Jo 
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i.e. mi = m 2 . Returning to (I4.18P we also have 

(4.24) f f {mi + m2){dxUi — dxU2f‘dxdt = 0, 

Jo Jo 

and so appealing to (j4.11l) and (|4.13p we have Ai = A 2 and Bi = i? 2 - Then it is straightforward 

to show that ui = U 2 by multiplying (j4.1jl by u = ui — ^2 and integrating by parts. Noting that 

Ux has an a priori bound, we have 

(4.25) -J u‘^{t,x)dx + — J J u^dxdt < C J J u^dxdt, 

and we conclude that n = 0 by Gronwall’s Lemma. □ 
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